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Abstract 

Using a prompt neutrino beam in which a Vr component was identified for 
the first time, the Vr magnetic moment was measured based on a search for 
an anomalous increase in the number of neutrino-electron interactions. One 
such event was observed when 2.3 were expected from background processes, 
giving an upper 90% confidence limit on of 3.9 x 10" ''/U^. 
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I. INTRODUCTION 



Magnetic moment measurements are an important tool for probing the fundamental 
structure of matter. Currently, precision measurements of the electron and muon magnetic 
moments are being used to probe the structure of the vacuum to the highest precision. 
A non-zero magnetic moment for any neutrino would be a clear and unambiguous signal 
for physics beyond the standard model. We have used the data collected for the Fermilab 
experiment E872 (DONUT) to perform a search for anomalous electromagnetic interactions 
of the tau- neutrino that would be a signature for a magnetic moment. This paper presents 
the first measurement of the Vj. magnetic moment using a neutrino beam in which the v^s 
have been positively identified. 

The direct limits on and magnetic moments are /i^^ < 1.8 x 10~^°/xb [|1|] and 
/ij,^ < 7.4 X 10~^°/iB and efforts are currently being made to extend these limits by 
another order of magnitude PQ. Because no direct measurements of the v-r component of a 
neutrino beam has been made until now, there has been no comparable measurement of its 
magnetic moment. In the CERN experiment WA66 0, a limit of < 5.4 x 10^ ''/i^ was 
based on an assumed flux of i^t-'s in the neutrino beam 0. Elsewhere, indirect limits were 
derived from the duration of the supernova explosion SN1987A, giving a limit of 10~^^^b 
for all neutrino flavors under the assumption that all neutrino flavors are equally produced 
in the supernova explosion. 

More recently it has been argued [0 that the oscillation of atmospheric z/^ to Ur would 
yield a limit of < 1.3 x 10"^//^. However, it has been pointed out that this value should 
be interpreted as a limit on the initial state i/^ instead 0. 

It used to be conventional to associate the properties of the neutrinos with their weak 
eigenstates z/g, t'^, and Ur- The discovery of neutrino oscillations P,p!0| implies that the 
electromagnetic properties of neutrinos should be associated with their mass eigenstates as 
opposed to their weak eigenstates . Since the parameters of the mixing matrix are as yet 
undetermined, and there remains the possibility that neutrino mixing is further complicated 
by an additional sterile neutrino |TT|, the composition of any neutrino beam in terms of 



the mass eigenstates can at this stage not be determined. Until the oscillation scenario is 
fully understood, the electromagnetic properties should therefore be characterized by the 
initial neutrino flavor. This description is independent of any oscillation assumptions and 
allows for an extraction of the mass eigenstate magnetic moments when all of the mixing 
parameters are known ||^. 

In our measurement of the tau-neutrino magnetic moment we searched for an anoma- 
lous increase in the elastic neutrino-electron cross-section above the value predicted by the 
Standard Model. In the Standard Model neutrinos interact with electrons through Z" ex- 
change, and a magnetic moment fi^, adds an extra component due to photon exchange. The 
cross-section for a neutrino interacting via its magnetic moment with an electron is given in 
the high-energy limit by 



dTe /i| m2 VTe E, 

where is the energy of the scattered electron in the laboratory frame [|1^. The total 
magnetic moment scattering cross section is obtained by integrating over Tg. The lower 
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integration limit is given by the experimental threshold for low-energy electrons (0.1 GeV in 
this analysis). Since the neutrino undergoes a spin- flip when a photon is exchanged, there is 
no interference with the Standard Model process, and the total neutrino-electron scattering 
cross-section is just given by the sum of the two contributions. 

Kinematic constraints limit the angle between the incoming neutrino and the scat- 



tered electron in the laboratory frame to be 

Bl_^ < 2mjE, (1.2) 

and for electron energies in excess of about 1 GeV, 6',^_e is less than 30 mrad. This angular 
constraint can be used as a clear signal to select neutrino-electron scattering events from 
the background of i/g-nucleon charged-current events in which the electron is produced at a 
much larger angle. 



II. THE APPARATUS 

The Fermilab experiment E872 (DONUT) took data in 1997 and z/t--N charged current 
interactions were observed in nuclear emulsion |T^. Having established the existence of z/^- 
in the u beam and measuring the ratio of z/g, z/^ and using charged current interactions, 
we searched for single electrons arising from u interactions, the signature of a z/ magnetic 
moment. The apparatus is described in detail elsewhere |T^; only the components central 
to this analysis will be discussed here. 

The experiment consisted of three essential parts: a Ur enriched neutrino beam produced 
in a beam-dump target, a neutrino target, and a spectrometer with electron and muon 
identification. A "prompt" neutrino beam was produced by a beam of 800 GeV protons 
incident on a tungsten target. The target length and material were chosen so that most 
of the long-lived secondary particles would interact in the target before decaying, while the 
short-lived particles would decay before interacting. Hence the main contribution to the 
high-energy neutrino flux came from the decay of the short-lived D-mesons, with u^s and 
z/^'s being produced in their semi-leptonic decays, and the z/^'s primarily from the decay 
channel Dg —>■ tVt, t Ur + X. Downstream of the tungsten target there was 30 m 
of absorber material and two sweeping magnets to remove all but the neutrinos from the 
beam. /^From our flux calculations using the LEPTO |]16| particle generator and a GEANT 
[ 171 simulation, we estimate that in the sample of located events,the neutrino interactions 



would be 47% z/^ charged current events, 27% z/g charged current events, 5% z/^- charged 
current events, and 21% neutral current events. This composition was confirmed by our 



measurements of charged current interactions [14 



Behind the absorber was the neutrino target region, shown in figure |l|, consisting of four 
steel and emulsion modules that were interleaved with planes of 0.5 mm diameter scintillating 
fibers. The entire region was surrounded by lead shielding, 20 mm thick upstream and 6 mm 
thick downstream. The total target mass, including the upstream lead, was 554 kg. In front 
of the lead shielding was a wall of scintillator counters used in the trigger to veto interactions 
upstream of the target region. 

The spectrometer consisted of the scintillating fiber tracker in the target region, drift 
chambers upstream and downstream of a wide-aperture analysis magnet, and a lead glass 
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array with a muon identification system behind it. The data acquisition system was triggered 
by a coincidence between a plane of scintillator counters just downstream of the target region 
and at least one of two planes of scintillators inside the region, together with no signal in 
the upstream veto wall. 

In the search for anomalous u-e interactions, the steel and emulsion target modules 
interleaved with the fiber tracker served as a high-granularity sampling calorimeter that 
allowed the identification and location of electromagnetic showers. 

III. EVENT SELECTION 

In the six-month run of the experiment 4.0 x 10^ events were recorded for analysis and 
2.0 X 10^ of these had two or more tracks in the target region. Candidate neutrino-electron 
scattering events were selected from these events if the following requirements were met: 
there were no identified muons associated with the event; there was a signal in the trigger 
hodoscope directly downstream of the reconstructed vertex; the reconstructed tracks had an 
angle of less than 100 mrad with respect to the neutrino direction, and the reconstructed 
vertex was in the fiducial volume (in a target module and no more than 0.24 m in the 
horizontal direction from the center of a module). 

These cuts removed events triggered by high-momentum muons that interacted in the 
material surrounding the target and its support system sending secondary particles into the 
target area. 

After these cuts were applied, events with a reconstructed vertex that was consistent 
with a neutrino interaction were selected in a visual scan, leaving 285 events in the data 
compared to 186 expected from a Monte Carlo simulation. The excess in the data was 
caused by the interaction of low-energy neutrons or photons in the most downstream target 
module. These events were removed by requiring a signal of at least 2 GeV in the lead glass 
calorimeter if the interaction vertex was in the most downstream target module. This cut 
reduced the number of events to 100 in the data and 95 in the Monte Carlo. 

The remaining z/-N scattering background events were removed if there were any tracks 
in the event not associated with an electron or if secondary particles from the nuclear 
breakup were identified upstream of the interaction vertex. Events were also rejected if 
any reconstructed track was identified in the muon ID system, if the ratio of the measured 
momentum of a track to the signal in the electromagnetic calorimeter was less than 0.5, 
if a track crossed more than two radiation lengths of target material without initiating a 
shower, or if the track left a sequence of large energy deposits in neighboring layers of the 
scintillating fiber tracker. Once these tests were applied 13 events remained. 

As a final test, events were selected if there was an identified electron at an angle of 
less than 30 mrad with respect to the neutrino direction and if there were no reconstructed 
large-angle tracks with angles in excess of 500 mrad. This left one event, shown in figure |1|, 
an interaction in the most downstream target module. 
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IV. DATA ANALYSIS 



The expected number of events from Standard Model processes was determined by a 
GEANT Monte Carlo simulation of the apparatus. The same program was used to find 
the selection efficiency for magnetic moment interactions. The accuracy of the simulation 
was tested by comparing data and simulation for two different sets of well- understood and 
easily identifiable events. The first set of control events consisted 200 z/^ charged-current 
interactions containing an identified muon. 

The second set consisted of 150 straight-through muons collected for calibration and 
alignment that scattered off an electron in the target region which in turn produced an 
electromagnetic shower. The electron energy spectrum of these "knock-on" electrons drops 
as l/Tg, while for magnetic-moment interactions it depends on l/Tg (see equation 
Though not identical, the similarity of the spectrum allows this control set to be used to 
quantify the effect of each selection cut on electromagnetic showers. 

As a test of our simulation each cut applied to the data was also apphed to both sets of 
control events. No significant difference between data and Monte Carlo was found with each 
cut removing the same fraction of events from both data and Monte Carlo for both control 
sets within statistical uncertainty. 

The Monte Carlo was also used to estimate the overall detection efficiency for electron 
events. As the event trigger required at least one hit in the most downstream of the trigger 
planes (T3), the shower from electrons produced in the neutrino interaction had to pass 
through two radiation lengths for each emulsion module. The acceptance was further limited 
by a 6 mm thick lead shield that was placed between the target region and T3. The overall 
trigger efficiency for events where neutrinos interact electromagnetically with electrons was 
(26.8 ± 1.4)%, ranging from 10% for interactions in the upstream lead wall to 50% for 
interactions in the most downstream module and, after the neutrino interaction selection 
cuts (15.2 ± 0.8)% of this type of events remained. By further requiring an identification 
of an electromagnetic shower in the event, the overall selection efficiency was reduced to 
(9.0 ±0.6)%. 



V. FLUX CALCULATION 



The fixed target D-meson production cross sections using a proton beam are crp{D 
(11.3±2.2)/i6/nucleon, ap{D°) = (28.0±2.5)/i6/nucleon, andcrp(D,) = (5.2±0.8)/i6/nucleon 
| |TE| . All cross sections were found to increase linearly with the atomic mass. For comparison, 
the previous CERN experiment used a value of ap{Ds) = 2.6/i6/nucleon at a beam energy 
of 400 GeV P], for which this calculation yields a value of 2.4/x6/nucleon. 

The differential cross section for D production was modeled by 

'^^'^ oc (1- |xi.|)"exp(-6p2) (5.1) 



dpfdxp 



with values of n = (7.4 ± 0.6) and b = (0.94 ± 0.06) GeV-^ 

The branching ratio for the decay Dg — rI7^ is (6.3 ± 0.5)%, based on recent measure- 
ments of the decay |T^. The resulting Ut fiux is (2.1 ± 0.3) x 10~^z/T-m~^POT^^. 
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VI. RESULTS AND DISCUSSION 



The single event surviving all the cuts occurred in the most downstream module and 
produced a narrow shower of particles with a reconstructed track at its center that had an 
angle of (10 ± 5) mrad with respect to the neutrino direction. The total recorded signal in 
the calorimeter was 20.0 GeV, of which 16.8 GeV were associated with the central part of 
the electromagnetic shower. For this energy, the electron angle should be less than 7 mrad 
according to Eq. (|1.2| ). The selected event could be a quasi-elastic z/g-N interaction with a 
relatively low-energy neutrino. 

The expected background rate due to neutrino-nucleon scattering was 2.3 events and 
other processes, such as neutral current neutrino-electron scattering, are expected to con- 
tribute less than 0.05 events. Since we observed one event when the expected background 
rate was 2.3 events, we conclude that we did not observe any signal events, and therefore the 
measured magnetic moment is zero. A statistical analysis based on the Feldman- Cousins 
method yields a 90% upper confidence limit of 2.3 events in the signal. 

To convert this limit in the number of signal events to a limit on the tau-neutrino 
magnetic moment, Eq. p..l| is integrated numerically, taking into account the calculated Vt 
energy spectrum. This yields 

= X 1.8 X 10-28m2. (6.1) 

In the experiment a total of 3.56 x 10^^ protons on target were recorded, and the average 
target mass during the run was 554 kg, corresponding to 1.7 x 10^^ target electrons, hence 
the number of expected events for a given magnetic moment, given by the product of the 
flux, cross section, and number of scattering centers, is: 

rievent. = ^ X 1.5 X 10^1 (6.2) 

This gives an upper limit for the magnetic moment of < 3.9 x 10^'' /i^- This should 
be compared to the experimental sensitivity of 4.9 x lO^^yU^, which is the limit that would 
be obtained had we observed as many events as predicted for the background. 

This analysis is flavor-blind and the limit applies in principal to the sum over all neutrino 
flavors. However, more stringent limits have been determined by other experiments for z/g 
and v^. Assuming a magnetic moment at the current limit, their contribution would be less 
than 10~^ events. The limit is therefore interpreted as a new upper limit on v^. 

While systematic uncertainties were not included in this analysis as they are not part of 
the Feldman- Cousins method, they would only change the result slightly, since the uncer- 
tainty due Poisson statistics completely dominates the estimate of the limit. Contributions 
to the systematic uncertainty come from the neutrino flux calculation (15%), the total pro- 
ton flux (15%), and the number of generated Monte Carlo events (5%). In addition to the 
statistical analysis based on the Feldman-Cousins method, a Bayesian analysis was per- 
formed including all these systematic uncertainties ||2^ using a flat prior distribution in the 
magnetic moment. This yields a 90% confidence limit of < 3.5 x 10~^/iB. We include 
this result for comparison. However, only the upper limit derived with the Feldman-Cousins 
method should be quoted. 
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VII. CONCLUSIONS 



The new upper limit for the tau-neutrino magnetic moment of 3.9 x 10''^ jjB is an im- 
provement over the previous hmit 0. Moreover, this is the first experiment to directly 
observe Uj. charged-current interactions and therefore be certain to have a Uj. component in 
the neutrino beam. 

The new limit is still three orders of magnitude above the limits for and t'^, and it 
dominates when extracting limits for the mass eigenstates from the current limits for the 
flavor eigenstates. Improving the limit on ^^^^ would require a u.,- beam that is comparable 
in total flux to previous z/g and z/^ beams. 
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FIG. 1. Target region view of the selected magnetic moment candidate event. The neutrino is 
incident from the left. 
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